Fermi surface reconstruction of superoxygenated La2Cu04 with ordered oxygen 

interstitials 
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Novel imaging methods show that the mobile dopants in optimum doped La2Cu04+y (LCO) get 
self-organized, instead of randomly distributed, to form an inhomogeneous network of nanoscale 
metallic puddles with ordered oxygen interstitials interspersed with oxygen depleted regions. These 
puddles are expected to be metallic, being far from half filling because of high dopant density, and 
to sustain superconductivity having a size in the range 5-20 nm. However the electronic structure of 
these heavily doped metallic puddles is not known. In fact the rigid-band model fails because of or- 
dering of dopants and supercell calculations are required to obtain the Fermi surface reconstruction. 
We have performed advanced band calculations for a large supercell Lai6Cu8032+jv where A'^=l 
or 2 oxygen interstitials form rows in the spacer LaieOie+jv layers intercalated between the Cu02 
layers as determined by scanning nano x-ray diffraction. The additional occupied states made by 
interstitial oxygen orbitals sit well below the Fermi level (-Ef) and lead to hole doping as expected. 
The unexpected results show that in the heavily doped puddles the altered Cu(3d)-0(2p) band 
hybridization at Ef induces a multiband electronic structure with the formation of multiple Fermi 
surface spots: a) small gaps appear in the folded Fermi surface, b) three mini-bands cross Ef with 
reduced Fermi energies of 60, 150, and 240 meV respectively, c) the density-of-states and band mass 
at Ef show substantial increases, and d) spin-polarized calculations show a moderate increase of 
antiferromagnetic spin fluctuations. All calculated features are favorable to enhance superconduc- 
tivity however the comparison with experimental methods probing the average electronic structure 
of cuprates will require the description of the electronics of a network of multigap superconducting 
puddles. 
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PACS numbers: 74.20.Pq,74.72.-h,74.25. Jb 



I. INTRODUCTION. 

The high-temperature superconducting cuprates 
(HTS) are heterostructures at atomic limit made of 
stacks of Cu02 atomic planes intercalated by spacer 
layers containing up to 25 percent of defects: oxygen 
interstitials, vacancies or atomic substitutions. The 
simplest cuprate La2Cu04 (LCO) is made of a single 
layer Cu02 intercalated by La202 rocksalt layers. The 
Fermi surface in La2Cu04 (LCO) is formed by a single 
band of the Cu02 layer that crosses the Fermi level 
[H . At half-filling the Fermi surface is destroyed by 
electronic correlations making a charge transfer gap 
between Cu3(f02p^ and Cu3d^'^02p^ many body 
electronic configurations. Insertion of dopants or defects 
in the spacer layers moves the chemical potential out 
off half-filling, forming a complex metallic phase. The 
rigid band model predicts a Fermi surface (FS) made 
of a large cylinder The breakdown of the rigid 

band model in the underdoped phase is shown by the 
appearance of i) Cu3(P02p^ states in the charge transfer 
gap 0; ii) a pseudo gap at the M point; iii) the Fermi 
surface made of short arcs Q and iv) an electron-like 
small FS pocket observed by quantum oscillations and 
high- field Hall effect [1]. In the underdoped regime the 
homogeneous metallic phase competes with electronic 
phase separation with the formation of a complex 



heterogeneous phase The charge-spin ordering 

and electronic phase separation are accompanied by 
fluctuations of the local structure that diverges from 
the average structure The defects inserted in 

the spacer layers to dope the copper oxide electronic 
structure are an additional source for the modulation 
of the heterogeneous electronic states near the Fermi 
level dl]. A particular attention has been addressed to 
oxygen interstitials that are mobile because of the misfit 
strain between the active and spacer lawrs inducing 
a contraction of the Cu-0 bonds [17|- Controlling 
their ordering by thermal treatments giant effects on 
superconductivity have been reported [l8l - [20| . Elec- 
tronic structure calculations show that the displacement 
of apical-oxygen, due to oxygen defects, changes the 
spin-phonon coupling[2l| and the electronic states near 
the Fermi level in cuprates [22|. The control of the 
functional electronic properties via dopants manipu- 
lation is of very high interest in manganites j23| and 
graphene [l^l- The focus has been recently addressed 
to the control of structural nanoscale phase separation 
into dopant rich puddles embedded in a dopant poor 
background in these complex materials. The dopant rich 
puddles of ordered dopants has been clearly observed 
also in iron chalcogenide superconductors and it has 
been shown that the two structural units have a quite 
different electronic structure with a superconducting or 
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a magnetic phase |25l - l27f . In the underdoped reghne of 
cuprates a structural phase separation between a dopant 
poor antifcrromagnetic phase and a dopant rich metalhc 
phase with doping close to 1/8 is clearly observed in 
La2Cu04+2; for 0< y <0.055 and YBazCugOe+y 
[29} . In the optimum and high doping regime of cuprates 
recent high resolution ARPES and STM experiments 
[30l - l32j show multiple electronic components with 
electronic phase separation at low temperature. A first 
pseudo-gap phase with doping around 1/8 competes 
with a highly doped metallic phase with about 1 /4 holes 
per Cu site is well described by band structure calcula- 
tions, confirming previous findings of phase separation 
giving distinct physics beyond the superconducting 
dome pl>,12] - Superoxygenated La2Cu04+y y>0.08 
provides an ideal case to study the phase separation in 
the optimum doping regime and the role of ordering 
of oxygen intcrstitials (Oi's) [33l - [36j . The O^'s in the 
rocksalt LaO layer are mobile because of the large tensile 
strain in the spacer layers. This is due to the lattice 
misfit between different layers [1^1 of Cu-0 bonds that 
is about 4 percent shorter than its equilibrium value 
of 1.97 A [131 • Therefore, by thermal treatments it is 
possible to induce their three dimensional ordering that 
has been shown to enhance the critical temperature to 
the highest value for La214 [s^. The variation of the 
critical temperature of the heterogeneous sample depend 
both on self organization of oxygen rich puddles and on 
the ordering of oxygen intcrstitials, as it has been shown 
using thermal treatments controlling of Tc [35j|. Using 
a novel experimental microscopy, scanning nano x-ray 
diffraction, probing both the real and the k-space the 
formation of a granular phase in La2Cu04+y has been 
discovered. Metallic heavily doped puddles, with ordered 
oxygen intcrstitials in the spacer layers, are embedded in 
a oxygen poor background (33 - [36| at optima doping. In 
the heavily doped grains the excess oxygen intcrstitials 
are ordered as shown in Fig. 1 forming stripes along the 
[1,1,0] direction. 

Although such compelling pieces of evidence for the 
importance of oxygen order there are no theoretical cal- 
culations on the electronic structure reconstruction near 
the Fermi level in the dopant rich puddles. In this work 
we present electronic structure results for the highly hole 
doped puddles of ordered oxygen interstitial in superox- 
ygenated La2Cu04. The method of calculation based on 
experimental structural information consider the super- 
cells of 0-rich LCO in sect. II. In sect. Ill we discuss 
the results of the calculations, and some ideas for future 
works are given together with the conclusions in sect. IV. 



II. METHOD OF CALCULATION. 

The calculations are made using the linear muffin-tin 
orbital (LMTO) method and the local spin-density ap- 
proximation (LSDA). Electronic correlation can be dis- 
regard because the oxygen rich puddles very far from 
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FIG. 1: (Color online) The ordering of oxygen interstitials 
(blue solid dots) in La214 in the (b,c) plane (upper panel) 
and in the (a,b) plane of the Fmmm structure as determined 
by x-ray diffraction using synchrotron radiation. The unit cell 
is Lai6Cug032+iv with N=2. 



half filling. The details of the methods have been pub- 
lished earlier [33]-[42- The elementary cell of La2Cu04 
(LCO) contains La sites at (0,0,±.721c), Cu at (0,0,0), 
planar O's at (0.5,0,0) and (0,.5,0) and apical O's at 
(0,0,±.358c), in units of the lattice constant ao=3.8 A, 
where c=1.24. In addition to the MT-spheres at the 
atomic sites we insert MT-spheres at positions (.5,0,±.5c) 
and (0,.5,±.5c) to account for the positions of empty 
spheres. The atomic sphere radii are 1.76 A (La), 1.26 
A (Cu), 1.18 A (planar O, and empty spheres), 
and 1.21 A (apical O), respectively. The z-projection 
of the corresponding irreducible Brillouin zone (BZ) is 
shown in Fig. 2^ by the limiting points F — A — M, 
whereas the IBZ for the doubled cell (La4Cu208) with 
ordered antifcrromagnetic (AFM) order would be lim- 
ited by F — A — Af" . Eight units of the elementary 
cell La2Cu04 are put together to form a long super- 
ccU, Lai6Cu8032, with the generating lattice vectors 
(1, — 1, 0), (4, 4, 0), and (.5,.5,c), so that its axis is ori- 
ented parallel to the (l,l,0)-direction, i.e. at 45 de- 
grees from the Cu-0 bond direction along (1,0,0). The 
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band calculations are made for this supercell containing 
72 sites totally, where one or two of the empty sites at 
the interstitial positions are occupied by excess oxygen 
ions. These supercells are chosen in order to represent 
fairly well the experimentally determined structures [sj] ■ 
A small lattice relaxation near oxygen interstitials is in- 
cluded in the first elementary segment of the supercell, 
so that apical-0 are pushed away about O.OSao in the 
x-y plane. The basis set goes up through i=2 for all 
sites. The z-projected IBZ corresponding to the super- 
cell is shown in|3^ by the hmits T - M3 - X3 - R. Self- 
consistency is made with 192 k-points and final results 
are based on 408 points in the IBZ. 



TABLE L Total and band decomposition of the DOS at Ef 
in units of {cell -eV)^^ , and average charges within each MT- 
sphere of the different atoms, in units of valence electrons per 
atom. Only 3 bands have a significant DOS at Ep, see the text 
for the numbering of the 3 bands in each case. The charges 
on Cu and especially La differ from site to site depending on 
the proximity to interstitial O. The O-charges are averaged 
planar and apical O and do not include the Oi's. 
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No interstitial sites are occupied with O in one set of 
the calculations (called LCO-0). The DOS of the LCO- 
calculations agree well with other published results on 
LCO fT]. The LCO-0 results serve as a reference for com- 
parison with the results with oxygen interstitials. For 
instance, the FS for (lightly doped) LCO forms a circle 
centered at the AI point in BZ of the elementary cell, 
but since the BZ corresponding to the supercell is folded 
and very flat along (1,1,0) it is necessary to identify the 
circle from several FS pieces in the folded zone. In two 
other sets of calculations we insert one (Lai6Cu8 032+i, 
LCO-1) or two (Lai6Cu8032+2, LCO-2) oxygen ions at 
interstitial positions. In LCO-2 the two O are put at 
equivalent positions in two adjacent AFM cells, so that 
two stripes of O^'s, separated by 1.41ao are formed along 
X (and z) [34. 



The excess O^'s sit at the interstitial interlayer posi- 
tions, above the oxygen ion in the Cu02 plane of the 
the orthorhombic unit cell, and form 3D ordered puddles 
below 350K, c.f. Fig. [TJ The insertion of Oi's is ex- 
pected to induce hole doping, because each new oxygen 
interstitial will bring 4 new bands well below Ep (one 
"s" and 3 "p"), but the oxygen has only 2 "s" and 4 "p" 
electrons. Therefore, simple arguments suggest that one 
Cu-0 band becomes unfilled, i.e. Ep has to go down rel- 
ative to the rest of the bands. However, this picture gets 
complicated by the fact that other atoms like La serve as 
charge reservoirs, lattice reconstructions are likely, and in 
addition the excess O positions are ordered in stripe-like 
patterns [3^ . 



It is often argued that correlation is too strong for 
having traditional bands in undoped cuprates when 
the Cu-d band is half- filled (43|. On the other hand, 
ARPES (angular-resolved photoemission spectroscopy) 
and ACAR detect FS's and bands that evolves with dop- 
ing and agree well with DFT (density-functional theory) 
calculations [H, 0,111]. Therefore, since we study doped 
cases here, it is not likely that correlation is more impor- 
tant than what is already included in our DFT calcula- 
tions. 
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FIG. 2: (Color online) Lower panel: The total DOS for 
Lai6Cu8032+jv with A^=0 and N=2. A broadening with a 
FWHM of ~ 0.01 eV is included. Middle panel: The dif- 
ference of the DOS for N=2 minus the DOS for iV=0. The 
maximum at -2 eV indicates the main energy location of the 
0(2p) impurity states of oxygen interstitial. Upper panel: 
The number of states, i.e. the integrated value of what is 
shown in the middle panel. 



III. RESULTS AND DISCUSSION. 

The nonmagnetic (NM) total DOS for Lai6Cu8032+Ar 
near the Fermi level are shown in Figs. [2J31 It can be 
seen that the DOS is increasing considerably near Ep 
when one or two Oi's are added in form of stripes. The 
DOS difference between the cases 7V=0 and N=2 in Fig. 
[^b shows that most of the Oi-p bands fall in the interval 
[-2.5,0] eV. Strong hybridization puts new states also on 
atoms around the Oi sites. The integrated differences (as 
the one shown in Fig. [^t) show that the electron count is 
about the same for all three cases at approximately -1.0 
eV below the respective Ep's. The DOS at Ep, N{Ep), 
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increases roughly three times when two O are added, see 
Table HI In the undoped sample the Fermi level Ep is 
well above the small peak in the DOS associated with 
the van-Hove singularity. In the doped cases the peak in 
the total DOS splits forming two or three peaks separated 
by a partial gap slightly below Ep, so that Ep is tuned 
exactly at the top of one the DOS peaks. The DOS 
functions from individual bands near Ep show similar 
shapes, but within a narrower energy range when the 
number of O^'s increases. The partial band DOS from 
bands 217-219 for N=2 are displayed in Fig. Eb; The 
shapes of the DOS from bands 211-213 for iV=0 look 
similar, but are ^3 times wider. Only three bands have 
a significant N{Ep), see also Table HI Of the two bands 
that contribute most to the total Ep, are the mid- and 
upper-band, 218 and 219 for N=2 (and 212-213 for iV=0, 
and 215-216 for 7V=1). These upper bands hybridize with 
Oi and make the local DOS on the added O very high, 
about twice as high as the local d-DOS on Cu. The local 
Cu-d DOS is largest on Cu-sites near the O^'s (typically 
twice the DOS on Cu far from O^). On the average there 
is a quite delocalized increase of the Cu-d DOS coming 
from the uppermost band when oxygen interstitials are 
added. This increase is smaller than of the total DOS, 
because of the large 0-p DOS, but it is most striking 
for the upper band, suggesting that this band is most 
important for any benefit of superconducting properties. 

The changes of the effective charges (see Table |l| show 
a weak hole doping on Cu when the number of O^ in- 
creases. A weak doping is also seen on the planar O 
charges in agreement with experiments showing only 
about one hole per Cu site induced by each oxygen in- 
terstitial at high doping. The charge on Oi is ~0.15 elec- 
trons less than within an apical O of the same volume. 
This difference makes the Is core levels to go up by about 
1.5 eV on O^ compared to the other O. 

The FS of undoped LCO consists of a rather 2- 
dimensional almost circular cylinder with origin at the 
M-point. The radius of the cylinder is quite large so 
that the FS englobes the M" -point in the center, see the 
schematic FS in Fig. |4^. The radius gets larger for in- 
creasing doping and the FS reaches the X- and Y-points 
at optimal doping when Ep coincides with the van-Hove 
singularity peak in the DOS. AFM order on Cu sites 
leads to a cell doubling and the BZ is folded along the 
X-Y line to obtain the irreducible part of the AFM-BZ 
within F-X-M" . One portion of the FS will be folded into 
the section " b" in Fig. |4^, while " a" is already within the 
AFM-BZ. Our band calculations are made for a 4 times 
extension of the AFM cell in the (1,1,0) direction, and 
the corresponding BZ is found from two more foldings of 
the AFM-BZ, the first at the line M2-X2 and once more 
at M3-X3. The new irreducible BZ is confined within F- 
M3-X3-R, and the portions "d", "c" and "e" constitute 
the down-folded FS circle of the elementary BZ into the 
BZ for the supercell. These 3 bands are found at Ep in 
the supercell calculations for LCO-0 and LCO-2 as can 
be seen in Figs. ^p-c. Fluctuations tend to open a gap at 
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FIG. 3: (Color online) Upper panel: The total DOS for dif- 
ferent bands near Ep in Lai6Cu8 032+jv for A'^=0,1 and 2, 
with a broadening of about 5 meV. Lower panel: The partial 
DOS of the 3 minibands crossing Ep in Lai6Cu8 32+2 and 
their sum (thin line) with broadening of about 3 meV. The 
Fermi energy of the upper miniband (219) is only 60 meV, to 
be compared with a superconducting gap of the order of 30 
meV. The Fermi energy of the other two minibands (218,217) 
are 150 meV and 240 meV. The Fermi energy of the corre- 
sponding bands (213-211 in the undoped case) for A'^^O is 
about 180, 300 and 950 meV, respectively. The gap between 
the 219 and 218 minibands, about 80 meV, is of the same or- 
der of magnitude as the so called pseudo gap in the anti-nodal 
point of the Fermi surface. 



the point "1" [42, |45|, |46| . It is also realized that portions 
"e" and parts of "c" will be gone at high doping. This 
is because the original M-centered FS circle will increase 
its radius for higher doping; the FS crossing that occurs 
on the M — X line for low doping will reach the AT-point 
at optimal doping, and ultimately cross the F — A line 
when the doping is increased further. The traces of the 
" e" and " c" portions will also be gone if there is a gap of 
the original FS at X, i.e. if the FS circle in the elemen- 
tary BZ breaks up into a banana-shaped feature centered 
at Af" . 

In Figs. HJd-c are displayed the FS pieces in the kz—Q 
planes of the irreducible BZ for N—0 and 2. For the 
calculations without interstitial O it is easy to recognize 
the FS-model in the real FS in the kz=Q plane. There 
are minor modifications in other kz planes because of 
a weak 3-D dispersion. The energy window is set nar- 
rower for the results with 2 interstitial O in order to 
avoid too many points for the more fiat bands in this 
case. Many of the same structures of the FS are seen 
in the plots for the results with interstitial O as for the 
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FIG. 4: (Color online) Upper panel (a): A schematic view 
of how a circular M-centered FS in the normal BZ (the ir- 
reducible part is limited by F-X-M) is folded into the BZ of 
the supercell given by the rectangle F-Ma-Xs-R. The M" is 
halfway between F and M and F-X-M" is the limit of the basic 
AFM BZ. Middle panel (b); the FS pieces (bands 211-213) in 
the fez=0 planes of the irreducible BZ (i.e. within F-M3-X3- 
R) for iV=0. Lower panel (c); the FS pieces (bands 217-219) 
in the fcz=0 planes of the irreducible BZ for N=2. 



results without them. But two differences draw our at- 
tention: First, there are many FS-points (this is more 
striking when the same energy window is used as for the 
undoped cell) to the left of the panel, close to F, when 
there are 0-interstitials, which suggests that mainly the 
FS piece "d" in Fig. |l]has gained a higher DOS. This 
agrees with the high DOS for the upper band (219 for 
N=2). In addition, section "d" is closer to the M3-X3- 
line, indicating hole doping as for an enlarged unfolded 
FS circle. Secondly, no FS is seen to the right (i.e at 



kp=0.5) for N—2, which probably means that a gap has 
opened near the X-point, or that the effective doping is 
such that the original unfolded FS already is beyond the 
X-point. A separate small hole pocket appears near ~ 
0.45. An effort to force a back folding of that pocket into 
the original zone would put that band close to X on the 
M"-X line. The possibility of a gap at X is interesting, 
since a pseudogap is often associated with enhanced Tc 
near optimal doping for conventionally doped LSCO, and 
a small pocket nearby sets up the condition for the Lif- 
shitz mechanism (47j . Such band features can to some 
extent be understood; A dispersive band in the [1,1,0]- 
direction breaks up by the new periodicity and the band 
becomes flat at the limit of the (shortened) BZ where 
even small gaps might appear. The DOS will be large if 
the band remains flat in other k-directions. It is a mat- 
ter of tuning the periodicity to the doping so that Ep 
coincides with a favorable structure in the DOS [41]. It 
has been demonstrated that the pseudogap depends on 
stripe-like potential modulations in the (l,0,0)-direction, 
as from phonon or spin waves [13, IH, , but here the 
stripes are oriented in the (l,l,0)-direction. These fea- 
tures are a bit different when only one oxygen interstitial 
is added. 

A high DOS is normally connected with strong spin 
fluctuations. It has been shown that for a certain corre- 
lation between the wave length of spin fluctuations and 
the doping, the pseudogap can open right at Ep, where 
also phonons can contribute |3^. Now the modu- 
lation is in the (l,l,0)-direction without optimization of 
the length of the super cell to the doping, but it is still 
of interest to see if AFM spin fluctuations can become 
stronger by interstitial O. This is seen from spin-polarized 
calculations where AFM fields are applied to Cu atoms 
to create a spin wave along the cell. Two pairs of Cu near 
and far from the O^'s have the applied fields {±5mRy), 
while other Cu at the nodes of the wave have no applied 
field. The local exchange enhancement (Stoner factor) 
and moments will be the measure of the strength toward 
spin fluctuations. 



TABLE II: Local Stoner enhancement, S, and within paren- 
theses the induced local moment (m in units of HB/atom) on 
4 Cu calculated for AFM fields of ±5mRy. The two first Cu 
are closest to the interstitial O, while Cu(3) and Cu(4) far 
from them are most disposed to show magnetism. Distortion 
in one layer of apical O and unequal distances to interstitial 
O in LCO-2 contribute to differences between Cu sites. 



Cu(l) 



Cu(2) 



Cu(3) 



Cu(4) 



(m) 



(m) 



(m) 



(m) 



LCO-0 1.46 (.061) 1.46 (-.060) 1.46 (-.060) 1.46 (.060) 
LCO-1 1.51 (.063) 1.42 (-.050) 1.62 (-.079) 1.63 (.081) 
LCO-2 1.45 (.060) 1.45 (-.059) 1.43 (-.055) 1.56 (.071) 



The results are summarized in Table HIl As seen, with 
extra O^'s the tendency is such that the local enhance- 
ments are larger on sites far from the oxygen interstitial. 
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on sites where the DOS is not highest. The DOS on Cu 
near oxygen interstitial is increased by direct hybridiza- 
tion with p-states on those 0-sites. The local 0-p DOS 
near Ep is large on O^ (larger than 1 eV~^) and nearby 
apical O (0.5-1 eV^^), which is in contrast to the case 
without interstitial O or far away from the 0—i where 
the DOS is small (less than 0.1 eV~^). This fact resists 
to an enhancement on Cu, because the Cu-d and 0-p 
bands are tightly bound and the latter are less adapted 
for magnetism. 

IV. CONCLUSION. 

Oxygen interstitials added in form of stripes in the 
spacer layers form three overlapping mini-bands crossing 
Ep making the DOS larger and Ep appears at the maxi- 
mum of one of the DOS peaks. Thus, the natural growth 
of Oi-stripes in the oxygen rich puddles seems to be an 
efficient way to generate potential modulations for hav- 
ing Ep on a DOS peak with high N{Ep), accompanied by 
enlarged electron-phonon coupling and Tc [4l| . There are 
many possibilities for spectroscopic verification of the lo- 
cal electronic structure near the interstitials. It would be 
very important to perform ARPES measurements using 
photon excitation beams focused on nanoscale spots to 
measure the electronic structure of the super-oxygenated 
puddles calculated here. It is clear that ARPES experi- 
ments using large photon beams give information on the 
spatially averaged structure probing at the same time 
both oxygen-rich and oxygen-poor domains. Therefore 
we propose that nano-focused photoemission experiments 
would be a powerful probe to detect band narrowing in 
oxygen ordered O^'s. The fact that we find an upward 



shift of the 0-ls core level of the order 1.5 eV on the sin- 
gle stripe of 0^ and nearby apical O (compared to the O 
far from the interstitials) suggests that the increase of lo- 
cal Oi-p DOS can be detected via the threshold energies 
of soft X-ray emission and/or absorption spectroscopy. 
The more than tenfold increase of the local p-DOS near 
Ep on interstitials is a large effect that should modify 
spectroscopic intensities as function of increased oxygen 
content. 

The changes of electronic structure on Cu sites are 
strong and moderately favorable to AFM spin fluctua- 
tions with possibilities for enhanced Xspin ■ Furthermore, 
there is a segmentation of the original FS with a super- 
conducting puddles, into smaller disconnected pockets, 
from which the shape resonance in the superconducting 
gaps near a Lifshitz transition emerges as an additional 
mechanism for high Tc JJ, 48] . We show that the ordered 
lattice with 3D oxygen interstitials in oxygen rich puddles 
in LCO provides a realization of a multi-band metallic 
puddles. All calculated features of the electronic struc- 
ture of oxygen rich nano-puddle are expected to be favor- 
able to enhance superconductivity however the compari- 
son with experimental methods probing the average elec- 
tronic structure of cuprates will require the study of net- 
works of superconducting multi-band puddles described 
in this paper. Finally theoretical works on HTS focusing 
on granular superconductors (49l . [soj . where the maxi- 
mum superconducting critical temperature occurs at a 
critical density for the superconductor-insulator transi- 
tion in a network of 2d percolation clusters [Ell could 
investigate networks of multi-band superconducting pud- 
dles for understanding many puzzles of high temperature 
superconductors . 
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